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Abstract 

I I In this work a physical modelling framework is presented, describing the intelligent, non- 

local, and anisotropic behaviour of pedestrians. Its phenomenological basics and constitutive 
elements are detailed, and a qualitative analysis is provided. Within this common framework, 
two first-order mathematical models, along with related numerical solution techniques, are 
derived. The models are oriented to specific real world applications; a one-dimensional model 

Sof crowd-structure interaction in footbridges and a two-dimensional model of pedestrian flow 
in an underground station with several obstacles and exits. The noticeable heterogeneity of 
the applications demonstrates the significance of the physical framework and its versatility in 
addressing different engineering problems. The results of the simulations point out the key 
^ role played by the physiological and psychological features of human perception on the overall 

^-H crowd dynamics. 
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m 1 Introduction 

cn 

In relatively recent years, crowd dynamics has stimulated the interest of scientists such as physi- 
cists, engineers, and applied mathematicians. This has been mainly caused by some media events, 
sometimes catastrophic, which have pointed out the necessity to complement traditional methods 
^ of scientific investigation, heavily relying on empirical observations, with prediction and simula- 

r '~j tion tools. For instance, in Civil Engineering it is well-known the case of the London Millennium 

rS Bridge, closed immediately after its opening because of excessive lateral oscillations developing 

when pedestrians crossed it. The bridge was reopened two years later after the installation of 
dampers to avoid uncomfortable vibrations. The intense research activity that has followed the 
event up to now (reviewed e.g., in [TJ[5j) found that those oscillations were due to unpredicted 
resonances triggered by synchronised walkers [5] and showed the importance to consider also the 
crowd-structure coupling when designing pedestrian facilities. Another significant example is the 
Jamarat Bridge in Mina (South Arabia), annual destination of Muslim pilgrims. The huge amount 
of people cramming the bridge produced in the Nineties hundreds of deaths for overcompression. 
In this case, mathematical models were successfully employed to simulate the flow of pedestrians 
and to investigate its relationships with the arrangement of the surrounding environment. This 
made it possible to highlight critical situations in which crowd disasters tended to occur and sug- 
gest effective countermeasures to improve the safety of the event f^. Aside from these serious 
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incidents, crowded places can be commonly experienced in daily life: shopping centres, railway 
or underground stations, airports, stadia, whose design may take advantage of virtual simulations 
to preliminary study and test different solutions, also in connection with safety and optimisation 
issues. This requires to set up suitable models of pedestrian behaviour, on one side accurate to 
catch the complex unsteady crowd dynamics and, on the other side, versatile and easy to handle 
to deal with real world applications. 

The modelling of pedestrian flows dates back at least to the Sixties, when experimental studies 
started to be flanked by ad-hoc models aiming at reproducing specific aspects of pedestrian be- 
haviour, such as queuing and route choice, see e.g., [5 and the review [S]. The first model suited to 
reproduce spatio-temporal patterns of pedestrian movement is probably [7], proposed in the first 
half of the Seventies. In recent years, many mathematical models have been proposed in the litera- 
ture to describe crowd dynamics, at both the microscopic ^ ^ iini HH lEl [13 and the macroscopic 
[m [m El HH US] scale. The physiological and psychological features of human perception are 
expected to play a key role in pedestrian behaviour. In particular, the visual perception implies 
that pedestrians acquire and make synthesis of non-local information about their surroundings, 
as widely investigated in the medical [20] and cognitive psychology [21] [22] fields. The non-local 
pedestrian perception suggests the need for developing non-local crowd models. Some features of 
visual perception (e.g., the visual field in the horizontal plane) have been considered in developing 
visually guided agent-based (microscopic) models of pedestrian behaviour in transportation engi- 
neering (see e.g., |13] and referenced papers) and in architecture at both building and urban scales 
(see e.g., [23]). A hint toward the use of space non- locality in first-order (i.e., mass-conservation- 
equation-based) macroscopic models can be found in [25] for vehicular traffic applications, in order 
to allow for instability phenomena characteristic of high density flows. Concerning pedestrians, 
some proposals to account for their anisotropic and non-local behaviour at the macroscopic scale 
are contained in recently proposed models [HI [H] . 

In this paper a phenomenological analysis is performed to point out the intelligent, non-local, 
and anisotropic features of pedestrian behaviour in normal (no panic) conditions. On these bases, 
a physical modelling framework is proposed. Then, two first-order mathematical models are de- 
rived, oriented to specific real world applications: a one-dimensional (ID) model to simulate 
crowd-structure interaction in footbridges and a two-dimensional (2D) model to simulate the 
crowd flow in an underground station with complex plan. These models resort to non-linear mass 
conservation laws with non-local fluxes, hence they directly deduce pedestrians' velocity from phe- 
nomenological principles, and ultimately provide a macroscopic description of the spatio-temporal 
crowd dynamics. The paper is organised into three more sections that follow the conceptual ap- 
proach described above. In Sect. [2] the phenomenological analysis, the physical model, and its 
qualitative analysis are described. In Sect. |3| two different mathematical models and the related 
numerical solution techniques are derived to handle ID and 2D crowd dynamics for real world 
applications. Finally, in Sect. |4] conclusions are outlined and some further research perspectives 
are briefly sketched. 

2 Physical framework 
2.1 Phenomenological analysis 

The proposed physical modelling framework is developed on the basis of a phenomenological anal- 
ysis of crowd dynamics, specially focussed on the role played by pedestrian perception, evaluation, 
and reaction to the surrounding conditions. In the following, the retained features of pedestrian 
behaviour are summarised in three conceptual blocks. 
Remarks about intelligent behaviour: 

(fl) pedestrians are active agents, i.e., under normal conditions without panic onset, they share 
the same objective of walking with the maximum velocity towards a target T (e.g., doors, 
exit, displays), bypassing possible obstacles and avoiding the most crowded zones (see the 
concept oi personal space developed in human sciences and psychology, e.g. in [15]). These 
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strategies enable them to determine actively their walking direction and velocity, without 
being passively subject to the laws of inertia; 

(f2) pedestrians are intelligent agents, i.e., their mind evaluates, selects, and/or makes synthesis 
of what it perceives according to various psychological criteria (e.g., the level of anxiety [17] 
or the capacity to perform ensamble evaluations [H]). 

Remarks about non-local behaviour: 

(f3) under normal external and subjective conditions (e.g., the area where pedestrians walk, 
hereafter called walking area, is illuminated allowing visual perception), pedestrians do not 
perceive the real world locally in space, due to their ability to see up to a given extent around 
them. Let us call sensory region "the area required by pedestrians for perception, evaluation 
and reaction" , previously introduced by Fruin in his pioneering work |28j , and let us notate 
it Rs] 

(f4) pedestrians react after a time interval has elapsed from the perception time (see e.g., [201 [52]): 
from a general point of view, a reflex time delay ti , which does not involve a long perception 
time nor synthesis, is needed for reaction to the actual conditions, while a volitional time 
delay T2 is required to scan the visual field and to adopt walking strategies. It is expected 
that Ti <^ T2 T , T being the reference macroscopic time scale required to cross the walking 
area at the maximum walking speed vm] 

(f5) from the previous features, it comes that pedestrians in a given position at a given time 
react to the conditions perceived in front of them at a delayed time, i.e., in a non-local way 
in both space and time. 

Remarks about anisotropic behaviour: 

(f6) pedestrians are anisotropic agents, i.e., they are not equally affected by stimuli coming from 
all directions in space. Specifically, they distinguish between ahead and behind, in normal 
situations being essentially sensitive to what happens in a symmetric visual field focused on 
their direction of movement; 

(f7) walking pedestrians adapt the depth and width of the sensory region to their travel purpose 
and to their walking speed (e.g., a pedestrian walking for leisure purposes is expected to 
scan a wider field than a commuter attaining a train, and the faster a pedestrian walks the 
deeper the space required to evaluate and react is). 



2.2 Physical model 

On the basis of the above mentioned features, a physical modelling framework is proposed. Let Q 
be a bounded domain representing the walking area (Fig. [T^), and assume that the target T lies 
on the boundary dfl. The proposed model is expressed by a mass conservation equation, herein 
written for a control mass M: 

where t is time. 

In order to fully describe the dynamics of the system, the velocity field of the moving mass 
should be determined as a state variable or be a priori known. First-order models relate the 
velocity field to the mass itself. 

Let us introduce the pedestrian velocity field v in the point x G £7 at time t: 

v(x, i) = i)(x,t)e„(x,i), (2) 

where v is the magnitude of v (i.e., the walking speed) and e^, is a unit vector representing the 
walking direction. In order to account for the contribution of the strategies outlined in point fl, 
the walking direction is modelled as the superposition of two contributions: 
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(a) (b) 

Figure 1: (a) Walking area il with obstacles (x=pedestrian position, T=target); (b) Sensory 
region R^, desired direction e^, and interaction direction 



(ml) the desired direction, i.e., the direction that pedestrians would follow to reach their target in 
a clear domain (the so-called unconstrained conditions [H]). It is basically affected by the 
geometrical configuration of the boundaries of the walking area e.g., perimeter walls and/or 
inner obstacles, and it does not evolve in time. This direction is identified by the unit vector 
ed (Fig. 

(m2) the interaction direction, i.e., the direction that pedestrians would follow to avoid crowded 
zones. This direction is induced by the interactions among pedestrians and is identified by 
the unit vector (Fig. [ij^), which, unlike ej^, evolves in time. 

It is worth pointing out that, in general, both directions may be the result of a perception process. 
In particular, both the geometry of the walking area and the crowding of the environment might be 
discovered by pedestrians during walking. Nevertheless, in the following pedestrians are assumed 
to already know the entire walking area from previous experiences, so that the desired direction 
is a priori determined. The knowledge of crowding and the related interaction direction are instead 
obtained through a perception process characterised by a non-local and anisotropic behaviour. 
The resulting walking direction e^, is expressed as: 



\9ed + (l-e)e,\' 

where 9 G [0, 1] is a dimensionless parameter that weights pedestrian attitude to give priority to 
the walking area layout or to the crowd conditions. It is worth noticing that reversed pedestrian 
flow occurs if ^ e [0, 0.5]. Notice that in the one-dimensional setting Eq. ^ simply reduces to 
e^, = i, where i is the unit vector in the direction of the horizontal axis, because the dimensional 
constraint forces both and to be directed along i. 

The walking speed v is generally determined as a function of the crowd density through one of 
the constitutive laws proposed in literature (the so-called fundamental relations reviewed in [29) 
or recently proposed in [4} 115)). In the following, the Kladek formula introduced by Weidmann 
[50] is adopted in the revisited version proposed in [5T] : 



wm < 1 - exp 



-7 



1 1 

Pp PM 



(4) 



where vm is the maximum walking speed in unconstrained conditions (also called free speed); pM 
is the jam density; 7 is an exponent, obtained through a fitting of the data in [281 132] (Fig- [2]), 
that makes the relation sensitive to different travel purposes (L = leisure/shopping, C = com- 
muters/events, R = rush hour/business); pp — Ppi'x., t) is the perceived density that a pedestrian 
located in x feels within his/her sensory region. 

Both magnitude and direction of the walking velocity require the modelling of the non-local 
anisotropic behaviour of pedestrians (points f3 to f7). 
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The sensory region Rg = i?s(x, e^, 5, a) is modelled as a portion of the ball centred in x with 
radius (5 and angular span 2(5, which reproduce the characteristic depth and width of the visual 
field of pedestrians, respectively (Fig. [T]d, point f7). Hence, in a 2D domain the sensory region is 
defined as: 

i?s(x, Gd, (5, a) = 1^ e f7 : 1^ - x| < (5, ' ed(x) > cosaj , (5a) 

In the ID case this definition simplifies as 

i?^(x, 5) = [x, a; + (5] C [0, L], (5b) 

the interval = [0, L] being the f D domain. The visual field is oriented according to the desired 
direction and is symmetric with respect to the latter, because of pedestrian attitude to bear 
in mind their target (point f6). Within the perspective of a simplified model, let us assume that 
the depth 8 depends on the walking speed only, while the width la is assumed to be constant. In 
particular, the radius 8 is expressed as a function of the local delayed walking speed: 

5(x, t) = (5(«(x, t - ri)), (6) 

where t\ is the reflex time delay introduced in point f4. It is worth noticing that b should be more 
properly computed with reference to the walking speed in the point x', namely the position at 
time t — T\ of the agent who at time t is in x, rather than in x itself. This essentially amounts to 
track pedestrian trajectories under a Lagrangian point of view, whereas Eq. ([6]) tries to keep an 
Eulerian reference. Such an approximation is however acceptable in view of the usual smallness of 
the reflex time delay ti, which, from the mathematical point of view, implies a small error |x' — x| 
for a smooth velocity field v. 

The constitutive law linking 8 to the walking speed v is expressed as follows: 

m - — « + <5o, (7) 

where As = As(x) and (5o > correspond to the maximum and minimum visual depth, respec- 
tively, along the direction in the empty walking area (Fig. [s]). The expression of (5(w) is inspired 
to the distance d of two successive pedestrians, defined at the microscopic scale in |33j on the basis 
of experimental results (Fig. |4|. In this perspective, and are understood as the macroscopic 
counterparts of c?(um) and c?0: respectively. In particular, 8q is assumed to be a constant equalling 
a given fraction of the reference length L. 
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Figure 3: Examples of As in 2D and ID spatial domains 
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Figure 4: d — u law proposed in [3 3) 



The visual field width a varies in the range [0, 90°]: in the following, the value a = 85° is 
adopted according to [23] . It is worth stressing that the visual focus is higher in the neighbourhood 
of a = and decreases for a approaching a |21j . 

The crowd density pp = /Op(x, t) and, in the 2D setting, also the interaction direction come 
from an intelligent evaluation process (points f2 and f6), which implies some weighting of the true 
density p. Different evaluation processes can be modelled, in order to describe the main features 
of the intelligent human behaviour within a crowd. In the sequel, four weighting strategies si, . . . , 
s4 will be presented, leading to different definitions of Pp and in terms of p in i?^ . All of them 
can be formulated following a common conceptual line: 

(i) Localisation. A point Xp G is localised by means of a functional relationship in space 
involving the true density p: 

where 'I' is some localisation operator in Rg. The point Xp has to be understood as the 
location where crowding is influential according to the specific perception criterion expressed 
by 

(ii) Perception of the density. In the strategies si and s2, the point Xp is used to evaluate the 
density Pp perceived in x as 

Pp(x, t) = p(xp, t), 

meaning that pedestrians are mainly affected by the crowd in the point that has drawn their 
attention. In the strategies s3 and s4, the perceived density is the result of an ensemble 
evaluation of the crowd distribution over a part of or in the whole sensory region. 

(iii) Route choice. Once pp is known, the walking speed v is computed from the fundamental rela- 
tion Q. Furthermore, in the two-dimensional case, the interaction direction is recovered 
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from Xp as 

the minus sign being due to that interactions lead pedestrians to steer clear of the crowd 
perceived in Xp. 

Four possible definitions of the localisation operator 5*, amounting to as many behavioural 
strategies from pedestrians, are now discussed, pointing out different psychological attitudes. 

(si) Pedestrians evaluate the perceived density at the intersection between the far boundary of 
the sensory region and the desired direction, independently of the density distribution: 

Xp = X + (5erf(x), (9a) 

which in the ID case reduces to fixing 

Xp — X + S. (9b) 

In this strategy, pedestrians are expected to be strongly determined to reach an attracting 
goal (e.g., pedestrians leaving a room with one exit), so that they are slightly interested in 
the surrounding conditions. Hence, an a priori localisation process is adopted, focussing on 
the far boundary of the sensory region along the desired direction. 

(s2) Pedestrian attention is drawn by the point in the sensory region where the maximum value 
of p is attained (in case of multiple points with the same maximal density, the closest one 
to X is chosen): 

= = argmaxpd, i). (10) 
ieRs 

An analogous expression holds in the ID case, up to maximising p{^, t) for ^ G [x, x + (5]. 

This strategy is intended to account for pedestrian anxiety about high crowd density (e.g., 
pedestrians leaving a room with multiple exits) , hence an a posteriori localisation is adopted, 
based on the observation of the environmental conditions. However, a pointwise criterion is 
used because pedestrian anxiety does not allow for any kind of ensemble evaluation of the 
crowd conditions. 



(s3) The point Xp is still given by Eq. (10 1, but the perceived density is a weighted average of 



the true density in x and in Xp based on the distance rp = |xp — x|: 

Pp(x, t) = (1 - g(rp))p(x, t) + g{rp)p{xp, t), (11) 

where 

ffW = -4-P + l- 

Since < rp < 5, the function g ranges in [0.2, 1] and is monotonically decreasing with rp. 
In particular, the maximum g{0) = 1 reflects the fact that a pedestrian located in the point 
of maximal crowd density gives it the greatest importance, so as to avoid contact with other 



people. The form of Eq. (Ill in the ID case is straightforward. 

In this strategy, pedestrians are still mainly concerned with the highest crowd density, but 
they are able to give it due weight with respect to its relative position. Hence, the evaluation 
of the perceived density implies self-consciousness from pedestrians, who establish a relation 
between themselves and the environmental conditions (e.g., commuters attaining their final 
destination in rush hours). 
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Figure 5: The family of functions (12a I 



(s4) Pedestrians evaluate Xp as the centre of mass of the whole crowd distributed in Rs , and use 
it as an indicator of the location of the mean value of the crowd density. However, when 
scanning the angular width of the sensory region they focus mainly on the frontal zone (cf. 
also [H]), hence the detection of the centre of mass is weighted by the angular position —a < 
a < a in Rg, see Fig. [ijj. Weighting is carried out by a function G : [—a, a] — ?> [0, 1], which, 
coherently with the symmetry of Rs around e^;, is assumed to be even and monotonically 
decreasing from G{0) = 1 to G{a) to render the lateral visual fading. Among different 
classes of weighting functions complying with these criteria, the following is considered: 

G{a) = 1 - (^-^y , (12a) 

where the exponent /i > may vary based on the travel purpose. In general, it can be 
noticed that the higher /i the wider the angular span of the visual field, see Fig. [5j Because 
of the lack of experimental data, the linear trend (/i = 1) is specifically retained here. 

The perception point and the corresponding perceived density are evaluated as 

Xp = = — / ^P(|, t)G(a(^ - x)) dt Pp{^, = T^' (12b) 
"m-Fis J \Rs\ 

where rufj^, \Rs\ are the weighted mass of the crowd in Rs and the weighted area of Rs, 
respectively: 

mfl. - J p{i, t)G{a{^ - x)) dt \Rs\= J G{a{i - x)) d^. (12c) 
Rs Rs 



In Eqs. ( 12b I, ( 12c I, a — — x) is the angle that the vector ^ — x forms with the desired 
direction of movement ed(x). The ID counterparts of these formulas are obtained by setting 
G = 1, because formally a — 0. Then Xp is the actual centre of mass of the crowd distributed 
in [x, X + S], and Pp{x, t) the actual average of p therein. 

The strategy just outlined relates to curious pedestrians: ensemble evaluation and synthesis 
are performed to both localise and evaluate the perceived density (e.g., in leisure or shopping 
activity) . 

Pictorial representations of the above localisation strategies are provided in Fig. [6) in a ID 
domain and for a given spatial distribution of the actual crowd density in the sensory region. 

2.3 Qualitative analysis 

The qualitative analysis of the presented model develops in two steps. First, the diffusive effect 
induced by the non-local features of the model is analytically demonstrated by using a model 
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equation in the ID spatial domain. Second, the results obtained by using the four localisation 
strategies are numerically evaluated and discussed for a given true density distribution in both 
the ID and the 2D case. 



2.3.1 Analysis of the diffusive effect induced by the non-local equations 



In order to study the effects introduced by the non-local in space fundamental relation, a model 
equation in the ID spatial domain is analysed. Let us consider the mathematical model given by 
the ID mass conservation equation, closed by a generic non-local fundamental relation: 




= 



(13) 



Let us assume for simplicity that the space dislocation is evaluated by means of the simplest 
strategy among the presented ones (i.e., si), so that Pp{x, t) = p{x + S, t). Assuming that the 
extension of the sensory region is small compared to the characteristic spatial length (i.e., the 
domain size L), a first order Taylor expansion of p gives: 



so that 



dp 

p(x + S, t) ^ p{x, t) -t- — ix, t)5. 

ox 



dp 

v{pp) « v{p) +v\p) — S, 



(14) 



(15) 



where ' will stand henceforth for derivative with respect to p. With the approximation (15), and 
introducing the non-negative function v^p) = —pv'(p)d (notice that is a decreasing function of 
the density, Eq. (Ill)), the mass conservation equation can be expressed in local form: 



dp 
di 



d_ 
dx 



(pvip)) 



d_ 
dx 



Kp) 



dp 
dx 



0. 



(16) 



It is worth pointing out that the function v [m^s ^] plays the role of an equivalent kinematic 



viscosity induced by the non-locality parameter 5. Multiplying Eq. ( 16 1 by p and then integrating 
on the spatial domain SI — [0, L] gives 



dp 

dt 



p dx 



d_ 
dx 



{pv{p))pdx 



d_ 
dx 



dp 

dx 



pdx — 0, 



(17) 



whence, setting q(p) = pv{p), a further manipulation yields 



1 d 
2dt 



p dx 



l'ip)^pdx 



d_ 
dx 



dp 



(18) 
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Integrating by parts the second and the third term at the left-hand side, with periodic boundary 
conditions apphed to p, Eq. ( 18 ) finally implies 





hence the spatial norm of the density (i.e., its "energy" ) decreases in time as it typically happens 
in diffusion phenomena. In other words, the spatially non-local model does not only allow the 
pedestrian intelligent behaviour to be described but gives also the mass conservation equation 
formally a parabolic trend, which for instance prevents the formation of unrealistic shock waves. 

Analogous conclusions are expected to hold qualitatively also for the other localisation strate- 
gies and in the 2D case. 



2.3.2 Effects of the different localisation strategies on the perceived density 



The effects of the different localisation strategies described in Sect. |2.2| are now evaluated in both 
the ID and the 2D setting by applying them to a test case with a given density distribution. In 
the following, the non-dimensional form of the model is considered, meaning that length, density, 
and speed are rescaled with respect to the reference values L, and vm. The instantaneous 
crowd density in the non-dimensional 2D domain VL = [0, 1] x [0, 1] is taken to be 



p{xi, X2) = Po + Apexp 



{Xi - Xi^c)^ + ix2 - X2,cf 



(20) 



see Fig. [7^: the reference constant value po > is perturbed about the point Xc = (x\^c, 2:2. c) 
by a bell-shaped curve attaining the maximum value Ap > 0, so that the density p ranges on the 
whole from po to po + Ap. The parameter > defines the width of the bell. For the distribution 
in the ID case over the non-dimensional interval f2 = [0, 1] the cross-section at x^ = X2,c is used 
(continuous line in Fig. [9|. The coefhcients in Eq. ( [20| ) are set to po = 0.25, Ap — 0.3, £ — 1/35, 
xi^c — 0.4, X2,c — 0.5. Finally, in the simulations the fundamental relation Q is set for Europe 
and rush-hour traffic, that is, pM — 6 ped/m^, vm = 1.69 m/s, and 7 = 0.273p7\/. It is worth 
noticing that, given the previous value of £, the density bump is quite narrow, which reflects 
a strong non-equilibrium in the density distribution. In this situation, the ratio S/£ should be 
generally much larger than 1, meaning that the sensory region might fully encompass the bump. 
Significant differences among the various perception strategies are therefore expected. 

The perceived density pp resulting from each of the localisation strategies si, . . . , s4 is shown 
in Fig. [8] for the 2D case and in Fig. [9] for the ID case. Instead, in Fig. [TO] the localisation 
mechanism underlying strategy si is pictorially represented in the ID case as an example. The 
following considerations are in order. 

(i) In strategy si the perceived density is constantly equal to po throughout the domain in both 
the ID and the 2D case (Fig. [8^), which implies that pedestrians are not affected at all 
by the denser crowding about Xg. The reason for this outcome is that the sensory length 6 
is everywhere quite large, because in a clear environment the maximum visual depth A^ is 
not obstructed by intermediate obstacles. As a consequence, the perception point located at 



the periphery of the sensory region (Eqs. (9a), (9b)) never falls within the density bump if 



the latter is sufficiently narrow. This is particularly evident in the ID setting: pedestrians 
located upstream the bump (x < Xc) tend to perceive the density condition downstream the 



bump (Fig. 10) 



(ii) In strategy s2 pedestrians are sensitive to the highest density they perceive within their 
sensory region R^. In particular, as soon as the centre Xc of the bump falls in they react 
to the maximum density in the domain, namely po + Ap. As a result, in the 2D setting a 
large area of O features a perceived density pp equal to po -I- Ap, possibly also in points far 
from Xc (Fig. Isb) according to the local orientation of the desired direction of movement 
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Figure 7: (a) Filled contours and cross-sections of the instantaneous true density, (b) Desired 
direction e,^ for the same test case 
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Figure 8: Filled contours of the perceived density pp 



(Fig. [7|d). Conversely, in a smaller area immediately downstream the bump and in front 
of the target T the perceived density coincides with the true one, for the anisotropy of the 
sensory region makes pedestrians insensitive to what happens behind them. The analogy in 
the ID setting is that the perceived density constantly equals po + Ap for x < Xc, then it 
decreases to po for x > Xc following closely the profile of the true density. 

(iii) In strategy s3 pedestrians are still affected by the highest density perceived within the sensory 
region, but in this case they are able to weight this information on the basis of the relative 
position of the perception point. The resulting 2D perceived density pattern (Fig. [Sj:) is 
similar to that obtained with strategy s2, however in those points x whose sensory region 
encompasses Xc the maximum po + Ap is attenuated, as it is averaged with po according to 



the distance Vp = |xc — x|. Specifically, using Eq. ( 11 ), pp results of the order of po+g{rp)Ap, 
which does not exceed po + Ap because g{rp) < 1. The ID case highlights in particular that, 
for X < Xc, the perceived density gradually increases as the pedestrian position approaches 
the location of the peak of the true density. Conversely, for x > Xc the trend is completely 
analogous to that of strategy s2. 

(iv) In strategy s4 pedestrians operate a global synthesis of the mass of people distributed in 
the sensory region Rs, indeed their perceived density corresponds to an integral average of 
the true density within Rg. In the 2D case, this average also takes into account the lateral 
visual fading and the consequent relative enhancement of the frontal vision. The resulting 
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perceived density pattern (Fig. |8p, notice the change of colour scale with respect to the 
previous cases) is more homogeneous than that produced by strategies s2 and s3 (in the 
sense that the gap between the minimum and the maximum value of pp in is much smaller 
than in the other two cases), nevertheless it is not constant like that returned by strategy 
si (that is, the bump of density in Xc is somehow felt). Notice however that the maximum 
of pp in SI is definitely lower than that of p as a consequence of the average process. The 
ID case clarifies that pedestrians approaching the bump feel a much lower density than the 
peak value, whereas those located immediately downstream the bump react actually to a 
density lower than the local one, since they perceive the dispersion of the crowd in front. 

Comparing in particular the last three strategies, it can be noticed how the steep gradients of 
Pp featured by s2 are gradually smoothed in s3 and dramatically smoothed in s4, as a consequence 
of less and less pointwise perception criteria. In general, all strategies seem plausible and suitable 
to address various applications implying different crowd traffic conditions. 



3 Real world applications 

In this section, the physical framework previously described is used to study the fiow of a crowd 
in two specific real world situations: along a vibrating footbridge (ID domain) and in an under- 
ground station (2D domain) with obstacles and multiple exits. The most convenient modelling 
methods and numerical approaches are developed accordingly, in the frame of macroscopic first- 
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Figure 11: Scheme of the crowd-structure coupled model 



order models. Numerical simulations are performed with the non-dimensional model: the crowd 
density and speed are rescaled with respect to pM and vm, while lengths are rescaled with respect 
to the characteristic dimension of the walking area L. The velocity-density relation Q is adapted 
for the case of Asia and rush-hour traffic: pM ~ 7.7 ped/m^, vm = 1-48 m/s, and 7 — 0.273pM- 

3.1 Case study 1: ID model of crowd-structure interaction 

A ID macroscopic model is suitable to describe specific kinds of crowd events in some walking 
facilities, such as dense crowds walking along walkway, corridors, and footbridges. In fact, the 
flow along such facilities is mainly one-dimensional because of their line-like geometry (i.e., the 
length is one or two orders of magnitude larger than the width of the walkway) and because 
pedestrians usually share the main objective of reaching a common target (for instance during 
particular events, such as opening days or demonstrations, and/or because of the specific function 
and location of the facility with respect to transport facilities) . 

If the aim is to model the crowd flow along a lively structure, such as a slender footbridge, 
crowd-structure interaction should be properly modelled. In order to do so, the so-called par- 
titioned approach in its differential form is applied to the coupled crowd-structure system, 
according to the framework schematised in Fig. [TT] The coupled model has been developed in 
several studies, summarised in [21[3S]: its detailed description is beyond the scope of the present 
paper, and only some fundative concepts are introduced in the following to allow the reader to 
interpret the presented results. 

With respect to the physical framework previously introduced, some additional modelling is 
needed in order to describe the action exerted by the structure subsystem on the crowd subsystem. 
The following assumptions are introduced: 

(i) the motion of the structure, expressed by the envelope of its lateral acceleration z{x, t), 
reduces the walking speed; 

(ii) pedestrians are not affected by the motion of the structure if the acceleration is under the 
threshold of motion perception Zc {'zc ~ 0.1 m/s^ in the following); 

(iii) pedestrians react to the motion of the structure non-locally in time, i.e., with a time delay 
fi that is expected to be greater than the time interval between two succeeding steps (in the 
following fi = 1 s); 

(iv) after pedestrians have stopped at time tg because of excessive acceleration zm {'zm — 2.1 
m/s^), a stop-and-go time interval T2 (t2 = 5 s) should elapse before they start walking 
again. 

Notice that fi and T2 can be ascribed to the reflex and volitional time delays mentioned in Sect. 
[2j respectively. However, in the present context they are referred to pedestrian perception of the 
platform motion by means of their feet and legs, rather than by means of the visual perception. 

The mathematical model describing the ID crowd subsystem is given by the mass conservation 
equation, closed by a constitutive law linking the walking speed v to the perceived crowd density 



13 



dp_^ djpv) ^ 
dt dx 



(21) 



where the function w(pp) is as in Eq. Q and g{z) = g{z{x, t — fi)) G [0, 1] is a corrective function 
with the following qualitative trend: 



ZM - Zc 



\i Zc < z < Zm f\{t <ts y t>ts + 1-2) 
if I > ZM Ait<tsWts <t<ts+ fa). 



(22) 



Model (21 1 must be supplemented by boundary conditions at a; = and x — L. 



The crowd-to-structure action takes place in two ways. First, the crowd-structure total mass 
m is constantly updated by adding the pedestrian mass rric to the structural mass rrig. Second, 
the lateral force F{x, t) exerted by pedestrians is expressed as a function of both p and z. The 
interested reader can refer to [5S] for a description of the macroscopic force model. 

Finally, the structure subsystem is modelled as a non-linear three dimensional (3D) damped 
dynamical system, whose equation of motion can be written as: 



ds 
dt 



if [,s] =F(p,i), 



(23) 



where s — s(x, t) is the structural displacement; and are the damping and stiffness operators, 
respectively; F{p, z) is the applied lateral force introduced above. 



3.1.1 Numerical approach 



In order to solve the crowd model described by the scalar non-linear Eq. (21), suitable numerical 



schemes for time and space discretisation are required. The classical explicit Euler scheme is 
adopted for the advancement in time, while the spatial derivative is discretised using a finite- volume 
approximation with first-order monotonic upwind scheme, according to Godunov's theorem 
Therefore, the discretised mass conservation equation at the j-th cell becomes: 



At 



Ax 



(24) 



where Ax and At are the uniform grid size and the time step, respectively, and 9j±i/2 Siie the 
numerical fluxes at the cell interfaces. 

In spite of the classical numerical approach adopted, some attention must be paid in dealing 
with the non-local features of the proposed model. In particular, the above-mentioned Godunov's 
first-order upwind scheme j37) requires two main adjustments. 

First, let us recall the well known Courant-Friedrichs-Lewy (CFL) stability condition: 



sup— 1^,1 <1, 



(25) 



which involves the convection velocity Vc- Under the assumption that, for all fixed t > 0, the 
function p is a diffeomorphism, the expression of Vc in the non-local model is 



Vc{p, Pp) = v{pp)+pv{pp) 



dpp dx 
dx dp' 



(26) 



In order to define the fluxes (?j±i/2 appearing in Eq. (24), the convection velocity must be nu- 
merically evaluated at the cell interfaces. Adopting a linear interpolation scheme to compute the 
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local density Pj+1/2 and the non-local density from their respective values in the cells j 

and j + 1, the discrete convection velocity at the cell interface j + 1/2 reads 

Vc, 3+1/2 = «(Pp.j+i/2) + Pj+i/2v'{Pp, 0+1/2)^^, (27) 



where App = Pp,j+i/2~Pp.j a-^d ^P = Pj+i/2~ Pj- Two important features of Eq. (27) are worth 
a mention: 

• the sign of the ratio App/Ap affects the sign of the convection velocity and the selection of the 
upwind node in turn. In other words, the sign of v^. cannot be determined by simply referring 
to the so-called capacity density value in the speed-density constitutive law illustrated in Fig. 

• the value of the ratio App/Ap affects the CFL condition. In particular, Vc may tend to 
infinity for vanishing Ap, which would make it impossible to select a discrete time step At 



fulfilhng the CFL condition (251 



A threshold e > is set on |Ap|: if \Ap\ > e then Eq. (27) holds, otherwise |Ap| = e is imposed, 
because pedestrians are supposed to be unable to evaluate differences in local density smaller than 
e. In the following e — 10^^ is chosen and, in order to save computational resources, an adaptive 
time step is adopted. 

Second, the numerical flux at the interface requires a non-classical definition in order to pursue 
entropic weak solutions : 

f min j'JIp''. Pp) if Pj < Pj+i 

niax g(p', p^) if pj > p^+i, 
^.7?G[o, 1] 

where p'' = (1 — ri)pj + ?7Pj+i and = (1 — ?7)pp, j + rjpp^ j+i arc linear interpolations of the values 
of p, Pp at the sides of the interface j -|- 1/2. 

The Finite Element Method (FEM) is employed for the space discretisation of the 3D structural 
multi-degree-of-freedom model, while its advancement in time is obtained by means of the classic 
Newmark method. 

The two subsystems are characterised by non-matching grids in space thanks to the differential 
partitioning adopted, while they share the same discretisation in time 

3.1.2 Results 

The proposed model is applied to the simulation of the crowd dynamics along a real world cable- 
stayed footbridge located at Toda City, Japan. The footbridge connects a boat race stadium to a 
bus terminal: at the end of the races, the bridge is crossed by the spectators who leave the stadium 
to reach the terminal. During these events, the crowd density attains 1.5 ped/m^ [55]: in these 
circumstances, the bridge deck experiences excessive lateral vibrations due to crowd-structure 
interaction. 

The crowd flow is simulated under three different conditions (referred to as "set-up" in the 
following) of increasing complexity: (i) along the motionless deck (i.e., only crowd dynamics is 
studied); (ii) with an imposed deck lateral motion (i.e., only the structure-to-crowd action is 
accounted for); (iii) in the complete two-way crowd-structure interaction. These simulations are 
aimed at discussing the influence of the non-local pedestrian behaviour on the crowd dynamics, 
the effects of the lateral motion of the platform on the crowd flow, and the effects of the non-local 
pedestrian behaviour on the structural response of the footbridge. Only one localisation strategy, 
namely s3, is adopted in this application. 



The selected inlet boundary condition (b.c.) for the crowd density (Fig. 12 1) exhibits a steady- 
state regime, with p = 0.17 (corresponding to 1.3 ped/m^), comprised between two transient 
regimes marking the beginning and the end of the stadium evacuation. It should be noticed that 
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the density decreases in a smoother way than the initial increase. This is due to the assumption 
that the stadium evacuation abruptly starts at the end of the boat race with a sudden increase in 
the crowd density, while the complete evacuation is expected to be smoother in time. 

The footbridge structure is characterised by a two-span continuous steel box girder, a two- 
plane multistay cable system, and a reinforced concrete tower. The total bridge length L is about 
180 m and the road deck width is 5.25 m. The deck mass per square metre is 800 kg/m^ and the 
damping ratio is approximately 0.7%. The tower and the deck are modelled with elastic beam 
elements, while each cable is modelled using a single truss element (Fig. [12}d). The damping is 
modelled via Rayleigh stiffness proportional damping. More details about the structural model 
can be found in [55] . 

The envelope t) of the lateral acceleration of the deck within set-up (ii) is obtained as 
z{\ — exp(— /3t))0i(a;), where z = 0.25 m/s^ is the maximum value of the acceleration, (3 = 0.02 



and 01 (x) is the first lateral mode shape in Fig. 12 d. The time history in a; = 0.3 is represented 
in Fig. \\%. 




Figure [T3| summarises some results referred to the set-ups (i) (Figs. 13i-c) and (ii) (Figs. [I3p - 
f) in terms of the instantaneous crowd density and speed at three time instants corresponding to 
different regimes: the transient one, in which the crowd is gradually filling the deck span (called 
"filling gradient"), the "uniform crowd" regime, in which an equilibrium condition is expected 
along the motionless deck, and the final transient regime, in which the crowd density gradually 
decreases at the footbridge entrance ("vacating gradient" in the following). For the set-up (ii). 



the instantaneous lateral acceleration of the deck is also plotted (Figs. 13i-f). Notice that in the 
interval \x,x] the acceleration exceeds the threshold value Zc- 
The following considerations can be synthetically made. 



In the filling gradient regime, no effects of the non-local model arise on either p or v. Indeed, 
the adopted strategy s3 implies that the (discrete) perceived density equals the local one 



whenever the latter is monotonically decreasing. In the set-up (ii) (Fig. 131), the effects of 
the motion of the platform on the walking speed are clearly visible in the span interval in 
which the lateral acceleration is higher than Zc'. the higher the latter, the lower the walking 
speed. 

• In the uniform crowd regime, local and non-local models provide the same results in set-up (i) 



(Fig. 13 d), as expected in equilibrium conditions. The combined effects of platform motion 



and non-local perception are highlighted in set-up (ii) (Fig. 13 3). It is worth recalling that 
the walking speed is affected both by the local platform motion and, in the non-local model, 
by the perceived density. As a consequence, for x < x the two models predict different 
walking speed values in correspondence of the same density value. Indeed, in the non-local 
model the walking speed depends on the density perceived in a sensory region which includes 
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filling gradient (t = 2.49) uniform crowd (t = 7.47) vacating gradient (t=\\ .95) 
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Figure 13: Crowd dynamics: instantaneous density and speed for motionless (a-c) and moving 
(d-f) platform 




Figure 14: Crowd-structure interaction: deck lateral acceleration at a; = 0.3 with local (a) and 
non-local (b) model 



part of the interval [x, x], where the platform acceleration exceeds Zc- For analogous reasons, 
along the interval [x, x] the platform motion and the anisotropic (forward) perception justify 
the scatter in the maximum (minimum, resp.) value of p (v, resp.) and their skewed spatial 
distribution in the non-local model. In particular, the crossing of the two speed curves at 
X « 0.23 implies that there the local density is equal to the perceived one. 

In the vacating gradient regime, the non-local model avoids the shock wave in the density 
distribution, coherently with the results of the qualitative analysis (Sect, 
the platform motion (set-up (ii). Fig. 
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2.3). In addition. 



reduces the magnitude of the convection velocity 
in both models, so that the vacating wave front is shifted backwards with respect to set-up 

(i) (Fig. im). 



Finally, the structural responses predicted by the local and non-local crowd model in set-up 
(iii) are compared in Fig. 14 in terms of time histories of the lateral acceleration at the point 
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X — 0.3. On the one hand, the structural responses predicted by the two models share the 
same general trend in time: in particular, the peak lateral acceleration is reached during the 
uniform crowd regime. On the other hand, the acceleration peak value predicted by the non-local 
model is approximatively 1.4 times greater than the one predicted by the local model. A detailed 
investigation of the reasons for this difference would require a full analysis of the adopted force 
model, which is beyond the scope of the present paper. An example of the relationships among 
the structural response, the pedestrian force, and the crowd density and velocity can be found in 
[35j . In the present context, it is worth stressing that this difference is relevant in the engineering 
practice of footbridge design, because the lateral acceleration affects both pedestrian comfort and 
footbridge serviceability [T]. 



3.2 Case study 2: 2D model of pedestrian flow in built environments 

To approach the modelling of 2D pedestrian flows, a technique recently introduced in [121 155] . 
which slightly departs from the direct application of hyperbolic equations, is adopted. In fact, 
2D hyperbolic equations would involve analytical and numerical difficulties, such as the correct 
definition of the convection velocity in presence of non-local flux and/or in domains with obstacles, 
that are instead bypassed by using the formalism introduced in the above cited papers. 

The starting point is the mass conservation equation ([T]), which in Lagrangian form reads: 

^ p{^,t)di^O, (29) 



dt 

where 7(f , t) is the configuration at time t > of the material sub-domain £ C n. In more detail, 
the flow map 7 is such that x — 7(X, t) is the position occupied at time t > by the point 



initially located in X, therefore Eq. ( 29 ) expresses the conservation of the crowd mass contained 
in £ along the spatio-temporal patterns of the system. The flow map is related to the velocity 
V = v(x, t) by 

^(X,t).v(7(X,t),.) ^3^^ 
j(X, 0)=X. 

A major difference with respect to the ID case is that in the present 2D setting it is possible, 
indeed almost essential, to take into account obstacles scattered in the walking area. These make 
some portions of the walking surface inaccessible to pedestrians, hence from the modelling point of 
view they can be understood as internal holes to the domain Q. The desired direction must then 
be modelled coherently with these constraints, so as to point towards the target T while bypassing 



intermediate holes. The assumption that pedestrians a priori know the walking area (Sect. 2.2 1 
implies that depends on the geometry of the domain only. In other words, the desired direction 
is not affected by either the perception process or the dynamics of the interactions, hence it can 
be regarded to all purposes as a datum of the problem. A possibility, which allows to deal with 
quite complex domains in a relatively easy manner, is to obtain as the (normalised) gradient 
of a scalar potential u = u(x) : il — > K satisfying Laplace's equation: 

Au = in n, (31) 

e.(x) ^ (32) 
|Vu(x)| 

Boundary conditions are required along dQ, including the internal boundaries delimiting the obsta- 
cles, in order to identify the attracting destination as well as the perimeter walls that pedestrians 
cannot cross. A detailed analysis of these issues can be found in [TO] . 



3.2.1 Numerical approach 

In order to solve numerically the equations of the model it is necessary to design a suitable scheme 
for the approximation in time and space of the density p. In |39j a discretisation method of the 
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mass conservation equation in the form ( 29 1 is proposed and proved to guarantee properties of 



conservativeness, positiveness, and spatial accuracy. Before outlining the main underlying ideas, 



it is worth pointing out that this method is not concerned with the solution of Eqs. (31), (32), 
which instead can be performed a priori by whatever technique suitable for elliptic equations (e.g.. 
Finite Element Method). 

Let us introduce a time discretisation t„ = n/St, n G N, for a certain (possibly adaptive) time 
step At > 0, then let us describe the action of the flow map 7 during one time step (n — > rt + 1) 
by means of a one-step flow map 7„ : f2 — > in the following way: if x is the position occupied at 
time t„ by the point initially located in X (in other words, x — 7(X, t„)), then 

7„(x) =7(X, t„+i). (33) 

A natural approximation of the one-step flow map, coming from an explicit Euler discretisation 



in time of Eq. ( 30 ) , is 

7„(x) «x + v„(x)At (34) 
for v„(x) = v(x, t„), as if the velocity were constant during one time step. Likewise, an explicit 



Euler discretisation in time of Eq. ( 29 ) produces 



Jpn+md^= J Pn{^)dt (35) 
where Pn(x) denotes (an approximation of) the density p(x, t„) and E is any subset of the domain 

n. 

As far as the spatial discretisation is concerned, let us assume that: 
• the domain is a rectangle. 



• the domain is partitioned with an orthogonal uniform Cartesian grid {Eij} consisting of 
Ni X N2 elements in the horizontal and vertical direction, respectively, each Eij being a 
square of characteristic edge's length Ax > 0, 

• each obstacle can be described as the union of several grid cells, so as to be encompassed by 
the partition. 

The density p„ is approximated by a piecewise constant function taking the value p^^- > in the 
grid cell Eij. An analogous discretisation is applied also to the velocity v„, which in general is 
not known exactly for it depends on the unknown pn through the interaction component e^. The 
approximate value of v„ in Eij is denoted v"- G . This induces a further discretisation on the 



one-step flow map (34), that ultimately becomes a piecewise translation on the numerical grid 
(because the velocity is now constant in each cell). Plugging the approximate densities into Eq. 
( 35 ) and choosing E = Eij finally yields 



k=l 1=1 



where | • | denotes the area of its argument. 



Equation ( 36 ) provides an explicit numerical scheme for determining the new coefficients 
{Pij^^}i,j from the knowledge of the coefflcients {p^ jij at the previous time step. The pro- 
posed spatial discretisation is accurate with respect to the level of refinement of the numerical 
grid provided the parameters Ax, At are linked by the following CFL-like condition [55] : 

At 
Ax' 



max — |v^;.|<l. (37) 



Condition (37) compares well with Eq. (25 1, with the remarkable difference that it involves the 
values of the flux velocity rather than those of the convection velocity, so that the latter need not 
be determined. 
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Figure 15: (a) The map of the underground station used for numerical simulations, (b) The 
desired direction field e^. (c) Inlet b.c. on the crowd density 



3.2.2 Results 

In this section the 2D model is used to simulate a crowd leaving an underground station. The 
main goal is to investigate the effect of the various perception strategies in unsteady conditions 
and in built environments with obstacles. Figure |15^ shows the map of the underground station 
considered in the simulations, the grey blocks representing pillars and handrails delimiting the 
platform and some corridors, respectively. It is assumed that the left, right, and bottom edges of 
the walking area are perimeter walls, that pedestrians cannot cross but for the exits at the end of 
the corridors. 

Although it does not reproduce any particular station, the map features some basic charac- 
teristics quite common to many real stations. Pedestrians getting off the trains access the station 
from the platform, then they aim at the exits placed at the end of three corridors. The corre- 
sponding desired direction field is obtained accordingly by means of Eqs. (32), (31), see Fig. 
15 3. Specifically, the potential u is set to the maximum value {u — 1) at the end of the corridors. 
Conversely, the normal derivative of u, namely the normal component of e^, is prescribed at the 
boundaries of the pillars and of the perimeter walls in such a way that pedestrians do not tend 
to cross the edges while possibly sliding along them. The direction of movement e„ is determined 
according to Eq. ([3|, with the parameter 9 set to 6* = 0.7 to avoid reversed flow. 

The unsteady inlet b.c. on the crowd density is plotted in Fig. |15| ;, the maximum density 
value being p — 0.02 (corresponding to 0.15 ped/m^): such a time history is intended to represent 
the arrival of a train at the station and its subsequent restart. It is worth noting that the low 
value of p is chosen in order to avoid the congested regime (p < Pca, Fig. [2j3) in the platform 
and along the passages, and to allow a sound comparison of the strategy effects on the free fiow 
approaching the corridors. 

Let us comment on the crowd density spatio-temporal patterns resulting from the perception 
strategies si, , 



s4 (see Fig. 16 showing 20 contour levels equispaced in the range [0.2, 1]). 



• Strategy si is the most fluent: pedestrians do not hesitate to follow the desired direction 
because, the vectors e^, being parallel (Eqs. ([8|, (9a|), it invariably results e^, = 
(Eq. ([3|). Consequently, the distribution of the flow is essentially dictated by the domain 
geometry, chiefly by the arrangement of the obstacles in the walking area. Pedestrians getting 
through the right passage systematically head for the upper corridor, those getting through 
the middle passage head for the middle corridor, and finally those getting through the left 
passage head for the lower corridor. No significant swerves occur in front of the accesses of 
the corridors, therefore the total emptying time is the smallest one. 

• Strategies s2 and s3, quite similar to each other, are on the whole the less fluent. The crowd 
continuously disaggregates and aggregates because of the interactions, forming density lanes 
separated by low density areas. For instance, pedestrians getting through the right passage 
split off in two branches pointing towards the upper and the middle corridor, respectively, 
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Figure 16: Instantaneous crowd density fields (perception strategies column-wise, times row-wise) 



and leaving an almost empty space between them. Instead, pedestrians getting through the 
middle and the left passage ignore the upper corridor, preferring the middle and lower ones. 
This ultimately causes part of pedestrians coming from left and right passages to access the 
middle corridor. As a result, a higher crowding globally occurs, which slows down the total 
emptying of the station with respect to strategy si. For reference, the estimated emptying 
time is approximately 1.5 times the one needed in strategy si. 

• Strategy s4 is qualitatively intermediate between si and the pair s2, s3. Compared to 
the latter, the flow is more ordered and homogeneous: no sharp switches between high 
and low density zones are observed, because of the averaging implied by this perception 
strategy. However, the crowd does not simply follow the desired direction like in strategy 
si. Indeed, pedestrians getting through the right passage split up again between the upper 
and the middle corridors, while those getting through the middle and the left passage head 
preferentially for the middle and the lower corridors. Nevertheless, the homogeneity of the 
flow at the access of these two corridors causes a relatively small crowding, which ultimately 
allows for a quicker emptying of the station with respect to strategies s2, s3. In this case, 
the emptying time is about 1.2 times the one needed in strategy si. 

In order to synthetically analyse the topological differences in the flows of strategies si and s2, 
the integral pedestrian flux Qi across the entrance li of the i-th corridor is introduced: 



Qi{t) = J qix2, t) -1(1X2, 



li 

where q pv is the vector-valued pointwise flux. It is worth pointing out that each corridor 



width is = 1 m (see Fig. 15 1), thus the integral fluxes Qi are directly comparable. Figure 17 
shows the time history of the fluxes, their values being normalised with respect to the (constant) 
integral flux imposed at the inlet boundary (corresponding to 0.22 ped/s). Data coming from the 
numerical simulation have been processed by a low-pass filter in order to highlight the long-scale 
dynamics. 

As far as strategy si is concerned, the following considerations can be proposed. 
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Figure 17: Time evolution of the normalised integral fluxes Qi across the entrances of the three 
corridors: (a) strategy si, (b) strategy s2 



The maximal fluxes in the corridors are attained sequentially in time, according to the 
increasing distance of each corridor from the platform. This is mainly due to that pedestrians 
keep the desired direction in the area from passages to corridors, hence the time needed to 
access the corridors is essentially a matter of distance to walk. Maxima are of different 
magnitudes, decreasing from the upper to the lower corridor, because the right passage 



serves the largest platform area (see Fig. 15 d) and, in addition, the flows coming from the 



three passages do not mix. In synthesis, the corridor fluxes in strategy si are confirmed to 
be completely driven by the desired direction and the passage locations. 

• The flux across the lower corridor features more than one local maximum, while the other 
fluxes are monotonically decreasing after a single maximum. This is caused by the delayed 
arrival of pedestrians getting through the left passage and walking along the bottom wall. 
The latter has a blockage effect on their velocity field, resulting in a reduced speed and an 
increased local density. 

As far as strategy s2 is concerned, the fluxes in the three corridors are very different because 
of the crowd mixing induced by this behavioural strategy. The effect of such a mixing is visible in 
several features of the flux time history, that are outlined in the following. 

• Local maximal fluxes have comparable magnitudes in all corridors, by far lower than those 
attained in the upper and middle corridor under strategy si. Moreover, fluxes are signifi- 
cantly different from zero for a longer time interval than in si. In other words, the mixing 
is responsible for a lower and more balanced exploitation of each corridor during a longer 
time. 

• The fluxes in the upper and the middle corridors reach their respective first local maximum 
at the same time t w 1.3 with the same magnitude. Indeed the flow getting through the right 



passage splits between the two corridors (see Fig. 16 at t = 0.5 and f = 1) as a consequence 
of the modelled pedestrian attitude to observe environmental conditions. Because of the 
same attitude, the flow getting through the middle passage swerves from the middle to the 
lower corridor and may induce the simultaneous local maximum in the latter flux (see Fig. 



16 at t = 1.00). 



After merging in front of the corridors (see Fig. 16 at t = 1 and t — 1.5), the flows getting 
through the middle and the left passages head together for the lower corridor. Then the 
local density at the access of the latter increases, causing subsequently a new swerve of the 
flow from the lower to the middle corridor. This corresponds to the second local maximum 
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of the flux in the middle corridor (t « 2.2) with simultaneous local minimum of the flux in 
the lower corridor. 

On the basis of the results obtained from the simulations, strategy s2 seems the most appro- 
priate for modelling pedestrian behaviour in an underground station under low density conditions, 
indeed it induces a realistic redistribution of the flow among all of the available corridors. In 
particular, from a design point of view it could be argued that the middle corridor is the most 
sought after by pedestrians, hence it is expected to play a key role in the emptying process of the 
station. 

4 Conclusions 

In this paper a phenomenological framework for the modelling of crowd dynamics has been pro- 
posed, based on the concept of behavioural strategy developed by walkers. Three main elements, 
leading to the definition of a behavioural strategy, have been identified: (i) the will to reach some 
destination] (ii) a perceptional process, which determines how pedestrians evaluate and synthesise 
the information coming from the environment; (iii) a decisional process, which leads pedestrians to 
the final route choice based on the information elaborated in the previous two steps. In particular, 
the perceptional process introduces two crucial features of pedestrian behaviour, non-locality and 
anisotropy. 

These ideas have been translated into a multidimensional mathematical model able to describe 
the macroscopic spatio-temporal patterns of pedestrian movement. The model is based on the 
mass conservation equation supplemented by a direct modelling of the velocity of the walkers, 
which incorporates the above-mentioned decision-based dynamics. In particular, the speed of the 
crowd is related to the macroscopic density perceived by pedestrians via a fundamental relation 
motivated by experimental data, whereas the direction of movement stems from the interactions 
among the individuals. 

The model has been used to address two specific real world applications: (i) crowd-structure 
interaction due to coupled dynamics of pedestrian flow and vibrating footbridge; (ii) pedestrian 
movement in a two-dimensional built environment (specifically, an underground station) with 
several obstacles and exits. Ad hoc computational techniques have been developed, in order to 
conveniently face the difficulties posed by some non-standard terms of the equations (e.g., the non- 
local flux of the crowd) as well as by the coupling with the structure in ID or by the interaction 
with the obstacles in 2D. 

The results of the numerical simulations have shown that the behavioural strategy plays a 
fundamental role in determining the dynamics and the overall evolution of the system. This 
implies that a correct characterisation of the behavioural strategy is a major modelling task for 
obtaining useful qualitative and quantitative information from the model. It should be noticed that 
the strategy, and in particular the perceptional process, may change dramatically depending on 
the travel purpose and the emotional state of the crowd. In the paper a few alternative perception 
strategies have been proposed, ideally referring to resolute, anxious, and curious pedestrians. 
These may apply to different categories of walkers with different travel purposes, e.g., travellers 
catching a train, commuters in rush hours, shoppers. 

Remarkably, panic has been out of the scope of the present work. The proper modelling of 
panic is probably not simply a matter of perception strategy, but involves a change also in the 
route choice and in the importance given to the desired direction of movement with respect to the 
interactions with other individuals. In [14' it is suggested that pedestrians entering a panic state 
tend to follow chaotically other individuals, dropping any specific destination, therefore they are 
mostly attracted towards areas of high density rather than seeking the less congested paths. In 
addition, in [H I15| some proposals for modelling the transition from regular to panic conditions 
in the framework of first-order models using speed-density laws are provided. All of these ideas 
may be duly included in the behavioural framework presented in this paper, and they certainly 
deserve further study also from the modelling point of view. 
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